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Abstract. We present the first speckle interferometric observations of R CrB, the prototype of a class of peculiar 
stars which undergo irregular declines in their visible light curves. The observations were carried out with the 6 m 
telescope at the Special Astrophysical Observatory near maximum light {V = 7, 1996 Oct. 1) and at minimum light 
{V = 10.61, 1999 Sep. 28). A spatial resolution of 75 mas was achieved in the if-band. The dust shell around R CrB 
is partially resolved, and the visibility is approximately 0.8 at a spatial frequency of 10 cycles/ arcsec. The two- 
dimensional power spectra obtained at both epochs do not show any significant deviation from circular symmetry. 
The visibility function and spectral energy distribution obtained near maximum light can be simultaneously fitted 
with a model consisting of the central star and an optically thin dust shell with density proportional to . The 
inner boundary of the shell is found to be 82 i?* (19 mas) with a temperature of 920 K. However, this simple 
model fails to simultaneously reproduce the visibility and spectral energy distribution obtained at minimum light. 
We show that this discrepancy can be attributed to thermal emission from a newly formed dust cloud. 
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1. Introduction 

The R Coronae Borealis (RGB) stars are a class of un- 
usual objects characterized by sudden declines in their vis- 
ible light curves as deep as AV ^ 8. They are extremely 
hydrogen-deficient and also carbon-rich (e.g. Asplund et 
al. 200C and references therein). The RGB stars are 



thought to undergo the formation of dust clouds in ran- 
dom directions, and it is believed that a sudden decline 
takes place, only when a dust cloud forms in the line of 



sight (Loreta |1934 O'Keefe |1939D . A newly formed dust 
cloud is expected to be accelerated by radiation pressure 
and to move away, expanding and dispersing over months, 
as the object gradually returns to its maximum visual 
brightness. The mechanism of the dust cloud formation 
and its temporal evolution are, however, still poorly un- 
derstood. Especially, the location of dust formatio n is in 
dispu te: fa r from the star, ^ 20Rj, (e.g. Fadeyev 1986 , 
198S, Feast 1996 ), o r very close to the photosphere, ^ 2R^ 



troscopic observations as well as theoretical progress on 
dust formation suggest t hat th e latt er sce nario may be 
the case (see, e.g. Glayton 1996, Feast 1997 ), however, no 
definitive answer is yet available. 

The spectral energy distributions (SEDs) of RGB stars 
exhibit infrared emission peaks around 6^8 fim. The IR 
excess, which accounts for typically 30% of the total flux, 
is constantly present, regardless of the visual brightness 
of the central star. Therefore, the IR excess originates not 
from a single newly formed dust cloud, but mainly from 
a group of dispersed dust clouds with temperatures of ap- 
proximately 600 - 900 K. For example. Walker et al. (1996) 



fit the infrared (spectro)photometric data of R GrB with 
a 650 K blackbody. The study of IRAS observations by 
Gillctt et al. ( 1986| ) led to the detection of an additional, 
very extended "fossil" shell around R GrB, whose diame- 
ter is as large as 18', with a temperature of ~ 30 K. Walker 



(1994) found such fossil shells for at least four RGB stars. 



Payne-Gaposchkin 1963 ). Recent photometric and spec- 
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Apart from the detection of the fossil shells, most 
of the observational results on the circumstellar environ- 
ment around RGB stars were obtained by photometry 
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Fig. 1. Visual light curve of R CrB based on the AAVSO 
data. The epochs of the speckle interferometric observa- 
tions are shown by the arrows 

Table 1. Speckle interferometric observations. Ac/AA: 
central wavelength and FWHM bandwidth of the filters, 
Nt- number of speckle interferograms acquired for the tar- 
get, N-R '. number of speckle interferograms acquired for the 
reference stars, T: exposure time of each frame, S: seeing, 
and p: pixel size 





1996 Oct. 1 


1999 Sep. 28 


JD 


2450358.2 


2451450.2 


V (mag) 


7 


10.61 


Ac/AA(/im) 


2.191/0.411 


2.115/0.214 


Reference star 


HIP80322 


HIP77743 


Nt 


167 


990 


Nr 


198 


1044 


r(ms) 


150 


80 


S'(arcsec) 


1.4 


1.1 


p(mas) 


30.5 


26.4 


Field of view 


8'.' 9x8'.' 9 


5'.'lx5'.'l 



and spectroscopy. Recently, Clayton & Ayres ( 2001 ) have 
revealed extended C ii A1335 emission around two RCB 
stars, V854 Cen and RY Sgr, by long-slit spectroscopy. 
However, such direct information on the spatial distri- 
bution of material in the vicinity of the central star has 
been very rare up to now. In this paper, we present high- 
resolution speckle interferometry carried out for R CrB at 
maximum and minimum light. The properties of the warm 
dust shell will be derived by simultaneous fits of the ob- 
served visibilities and SEDs using power-law models. We 
will also discuss the possible indication of a newly formed 
hot dust cloud. 
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Fig. 2. a and b Two-dimensional power spectra of R CrB 
obtained on 1996 October 1 and on 1999 September 28, 
respectively, c and d Azimuthally averaged visibility func- 
tions derived from a and b, respectively 



2. Speckle interferometric observations 

The i^-band speckle interferometric observations were 
carried out with the 6 m telescope at the Special 
Astrophysical Observatory (SAO) in Russia, using our 
NICMOS-3 camera for the observation in 1996 and our 
HAWAII array speckle camera in 1999. The observations 
are described in Table 0. 

Fig. |l| shows the visual light curve of R CrB in the rel- 
evant period based on the compiled data of the American 
Association of Variable Star Observers (AAVSO). The 
data before 1996 were taken from AAVSO Monograph 4 
Supplement 1, while J. A. Mattei kindly provided unpub- 
lished data for the rest of the period. On 1996 October 
1, R CrB was on its final recovery from the 1995-96 deep 
minimum and its visual magnitude was around 7 (Mattei 
2000). Since V is approximately 6 at maximum light, the 
star was slightly obscured by AV = 1. On 1999 September 
28, the star was just between two sharp deep minima and 
had a visual magnitude of 10.61 (see Sect. ||). The star 
was heavily obscured by AV = 4.6. 

Fig. ^ shows the two-dimensional power spectra and 
azimuthally averaged visibility functions reconstructed 
with the speckle interferometry method (Labeyrie 1970). 
The error bars include systematic errors caused by see- 
ing variations as well as speckle noise errors. There is 
no significant deviation from circular symmetry in the 
two-dimensional power spectra observed at both epochs. 
The granular features seen in Fig. ||a are the speckle 
noise resulting from the small number of interferograms. 
In Fig. ||b, a slight elongation along the vertical axis of 
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the figure can be marginally recognized. However, this is 
within the error bars, and cannot be regarded as a sig- 
nificant feature. Clayton et al. (1997) suggest a bipolar 
geometry for R CrB based on the wavelength dependence 
of the position angle of polarization. Our speckle observa- 
tions seem to support a rather symmetric distribution of 
material. However, a disc or torus may lie very close to the 
star, and therefore, it is premature to rule out a bipolar 
geometry completely. 

The reconstructed visibility functions shown in Figs. |^ 
and ^ exhibit no difference larger than the error bars. 
The source probably consists of the central star (a point 
source) and an extended dust shell. Visibilities observed 
for such objects have a plateau at high spatial frequen- 
cies resulting from the point source. In the observations 
presented here, the angular size of the shell is small, and 
therefore, the plateau is not visible within a cut-off fre- 
quency of ^ 13 cycles/arcsec. The determination of the 
physical parameters of the shell by simultaneous fits to 
the observed visibilities and SEDs will be presented in 
Sect. H 

3. Photometric data 

UBVRIJHKLM photometry was carried out on 1999 
September 29, just one night after the speckle observation. 
The optical (UBVRI) data were obtained with the 1.25 m 
telescope at the Crimean Astrophysical Observatory, and 
the infrared {JHKLM) data with the 1.22 m telescope at 
the Crimean Laboratory of the Sternberg Astronomical 
Institute. 

No photometric data, except for the visual magni- 
tudes compiled in the AAVSO database, have been pub- 
lished around the date of our speckle observation in 
1996 October. Therefore, we use the data published by 
Shenavrin et al. ( 1979| ), who report a series of photomet- 
ric observations during the 1977 deep minimum. We use 
the UBVRJHK data on 1978 January 22 (JD2443530.6), 
when the visual magnitude was 7.03, almost the same 
as that on the date of our speckle observation in 1996 
October. For the L- and M-bands, we use photomet- 
ric data obtained on 1996 June 11 (JD2450246.4) by 
Shenavrin (unpublished observation) , about 4 months be- 
fore our speckle observation. It is known that the variation 
of the L- and M-band fluxes has no correlation with the 
visual light curve, and that they exhibit a semi-periodic 
variation of ~ 1.5 mag with a period of about 1260 days 
(e.g. Feast et al. 1997). In fact, the L and M magnitudes 
on 1996 June 11 are by 1.09 mag and 1.25 mag brighter 
than those on 1999 September 29, respectively. The L- and 
M-band fluxes measured 4 months before are expected to 
represent the values on the date of the speckle observation 
rather well. 

We also use IRAS observations at 12, 25, 60, and 
100 ^m. The very extended "fossil" shell gives rise to flux 
excess at 60 and 100 um, when evaluated with a large aper- 
ture (Gillett et al. 1986 ). In this paper, however, we are 
primarily interested in the warm dust shell which was not 



Optically thin dust shell 



Central Star r._ 



Cloud 




Obs. 

Fig. 3. Schematic view of the central star, an obscuring 
dust cloud, and an optically thin dust shell. The figure is 
not to scale 

resolved with IRAS, therefore, we adopt the point source 
processed fluxes for 60 and 100 /im. Long-term variations 
at wavelengths longer than 10 /xm are considered to be 
smaller than those in the L- and M-bands. Forrest et al. 



(1972) show that the amplitude at 11 /xm is about half 
as small as those in the L- and M- bands. The varia- 
tion at wavelengths longer than 11 /xm is estimated to be 
~ 0.75 mag at most. 

The interstellar extinction toward R CrB is E{B — 
V) = 0.05 (Asplund et al. 1997). The observed flux es are 
dereddened using the method of Savage & Mathis (1979) 
with Av = 3.1E{B-V). 



4. Description of tlie model 

The basics of the mode l used here are the same as adopted 
by Gillett et al. ( 1986 ). Fig. || illustrates the picture con- 
sidered in our modeling. It consists of the central star and 
an optically thin dust shell responsible for the constant 
IR excess. As mentioned in Sect. |2[ the star was obscured 
by AV = 4.6 (1999 September 28) and AV = 1.0 (1996 
October 1). Therefore, an additional obscuring dust cloud, 
which is assumed to be circular (projected onto the sky as 
seen from the star), is also considered, as shown in Fig 

Adopting Mboi and T^e to be -5.3 (G illett et al. p^86| 
and 6600 - 6900 K (Asplund et al.|2000D, respectively, the 



radius of the central star is approximately 70 Rq. We 
adopt 6750 K for the effective temperature, and radiation 
from the central star is approximately represented by a 
blackbody. T he dis tance to R CrB adopted here is 1.6 kpc 
(Gillett et al. 1986), yielding an angular radius of the cen- 
tral star of 0.23 mas. 



4 



K. Ohnaka et al.; DiflFraction-limited speckle interferometry of R CrB 



We consider a spherical shell with an inner radius rjn 
and outer radius rout • The shell is assumed to be optically 
thin not only in the infrared but also in the optical, since 
there is no evidence for strong circumstellar reddening for 
R CrB. Asplund et al. (1997) found that the observed op- 
tical flux at maximum light can be well reproduced by 
their line-blanketed model atmospheres without any cir- 
cumstellar extinction. In fact, the optical depth of the dust 
shell in the t^-band is ^ 0.3 in our models presented be- 
low. The assumption of spherical symmetry for the shell 
is justified, given the negative detection of deviation from 
circular symmetry in the observed two-dimensional power 
spectra. 

The grain number density is assumed to decrease in a 
power-law of radius. 



n{r) — Cr 



(1) 



where C is a constant which does not affect the shape of 
emergent SEDs, and therefore, should be adjusted to fit 
the flux observed on the earth. 

The temperature in an optically thin dust shell is de- 
termined by the thermal balance equation 



47rr-^ 



47ra^(5abs(a, v)Tr B^{T(r))diy, 



(2) 



where L^, is the luminosity of the central star at a given fre- 
quency, a is the radius of a grain, (5abs(a, ty) is the absorp- 
tion efficiency of a grain, and Bi,{T{r)) is the Planck func- 
tion. Amorphous carbon is the most probable candidate 
for th e circ umstellar dus t aro und RGB stars (e.g. Holm 
et al. |1982[ Hecht et al. [19841) . In Fig. |, we show tem- 



perature distributions predicted by the thermal balance 
equation, using the extinction of amorphous carbon ob- 
tained by Bussoletti et al. ( 1987] ) (AC2 sample). R ouleau 
& Martin ( |l99l[ ) (ACl sample), and Colangeh et al. ( [1995D 
(AGAR sample). We calculate Qabs from the complex re- 
fractive index derived by Rouleau & Martin (1991) in the 
Mie theory, assuming a single grain size, a = 0.01 /im, and 
using the code published by Bohren & Huffman ( 1983| ). 
This grain size is based on the result obtained by Hecht 
et al. (1984), who analyzed the UV spectra of R GrB and 
RY Sgr and concluded that the grain size is between 0.005 
and 0.06 /im. As Fig. ^ shows, the temperature distribu- 
tions agree with one another within ~ 50 K. 

The flux density observed on the earth from the opti- 
cally thin shell described above can be calculated by 



/s(A) = 



Bx{T{r)y-^ dr, 



(3) 



where and k\ is the mass and mass absorption coef- 
ficient of a grain, respectively, and D is the distance to 
the star. The constant C is adjusted so that the fiux pre- 
dicted from the models can reproduce the observed flux at 
25 /im, since its long-term variation mentioned in Sect. ^ 
is expected to be minimal, and at the same time, it does 
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Fig. 4. Temperature distributions derived from the thcr 
mal balance equation. Dotted line: Golangeli et al. (1995) 
dashed line: Bussoletti ct al. (1987), solid line: Rouleau & 



Martin (1991) with a = 0.01 /im 



not exhibit the contribution of the fossil shell. The varia- 
tion of the 25 /tm flux is not available in the literature, but 
it is estimated to be smaller than the 0.75 mag at 11 /im. 
Tout is set to be ~ 1.5 x lO'' Ri,, which gives a good match 
to the far-infrared part of the observed SEDs, as we will 
show below. 

It should be stressed here that the real circumstellar 
environment around R GrB is most likely much more com- 
plex than depicted by the above spherical shell model. 
Since dust ejection in RGB stars presumably occurs in 
clouds, not in a spherical shell, it is likely that the real 
distribution of material is clumpy or patchy without any 
clear inner boundary. At large distances, however, it is 
still plausible that a spherical shell model roughly repre- 
sents the distribution of material, as long as dust clouds 
are ejected randomly and frequently. There may be some 
newly formed clouds inside the inner boundary defined in 
the above spherical model, if the star ejects dust clouds 
frequently, for example, every pulsational cycle (40 - 50 
days). Without detailed knowledge about the dispersal 
processes of clouds, however, it is beyond the scope of this 
paper to construct a more detailed model. In the frame- 
work of our models, we consider only one newly formed 
dust cloud, as Fig. ^ illustrates. 

The observed SEDs (see Sect. % demonstrate that the 
contribution of the central star is not negligible in the 
near-infrared. The fiux from the central star was atten- 
uated by a dust cloud in front of the star by /S.V ~ 1 
on 1996 October 1 and /W - 4.6 on 1999 September 28. 
In order to estimate the attenuated flux from the star in 
the near-infrared, the extinction curve of the dust cloud is 
empirically derived from photometry in the optical, where 
emission from the dust shell is negligible. The effect of the 
obscuration due to a dust cloud is expressed as 



/4A) = /r'^(A)exp(-T^'(A)), 



(4) 



where /^(A) and f^'^^{X) denote ffuxes observed at any 
given time and at maximum light, respectively, and t'^'(A) 
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Fig. 5. Extinction curves of a dust cloud and its tem- 
poral variation. The open symbols represent the extinc- 
tion curves derived from the observations covering the 
1983 minimum obtained by Goncharova ( 1992| ). Open cir- 
cles: 1983 Oct. 17 (JD2445625.21), open squares: 1984 
Jan. 16 (JD2445715.60), and open triangles: 1984 May 22 
(JD2445843.43). The filled squares represent the extinc- 
tion derived for 1999 Sep. 29. The filled triangles repre- 
sent the extinction derived from the data obtained on 1978 



Jan. 22 by Shenavrin et al. (1979). The visual magnitude 
at each epoch is also given. The extinction characterized 
by 1/A is plotted with the dotted line for reference 



is the optical depth of the dust cloud in front of the star. 
Note that the effect of interstellar extinction, though small 
for R CrB, cancels out, because it affects both /*(A) and 
jmaxj-^-j ^jjg same amount. 

Fig. 1^ shows the temporal change of the optical depth 
of a dust cloud during the 1983 minimum in the U, B, 
V, R, and / bands. The optical depth in each band was 
derived using the photometric data in the decline and at 
maximum light obtained by Goncharova (1992). We also 
plot the extinction curves derived for the photometric data 
discussed in Sect. |[ The figure illustrates that the ex- 
tinction is almost independent of the wavelength at the 
very bottom of the deep minimum {V = 14.19). As the 
star starts its fina l reco very, the extinction curve starts to 
steepen. Pugach ( 1984 ) proposed that the neutral extinc- 
tion observed during the initial drop to minimum light 
can be explained by an optically very thick dust cloud 
whose coverage over the stel lar di sk varies with time. On 
the other hand, Hecht et al. (1984) analyzed the change of 



the ratio of total to selective extinction {Av/E{B — V)) 
during decline events for R CrB and another RGB star, RY 
Sgr, and suggested that rather large glassy carbon parti- 
cles with radii from 0.075 to 0.15 /^m might be responsible 
for the nearly neutral extinction seen at the very bottom 
of minimum light. These sizes are significantly larger than 
the 0.01 /im we adopt for a single grain size, but it should 
be noted that we adopt 0.01 /im as the grain size in the 
optically thin dust shell, while the large grains mentioned 
above are claimed to be present in a newly formed op- 
tically thick dust cloud. For the temporal change of the 
exti nction at the rise phase from minimum light, Hecht et 
al. (1984) proposed shat tering of large grains as a possi- 
ble mechanism. Efimov (1990) also proposed that the ob- 
served temporal variation of colors and polarization can 
be accounted for by the change of grain sizes, but assum- 
ing graphite grains instead of amorphous carbon. In any 
case, the spectral index of the dust cloud approaches 1, as 
the cloud disperses and becomes part of the optically thin 
dust shell. Therefore, the dust properties in the optically 
thin shell may presumably remain constant. 

The lack of detailed knowledge about the properties 
of dust grains formed in R CrB and of their tempo- 
ral variation forces us to adopt an empirical extinction 
law. As Fig. II shows, the extinction due to a dust cloud 
can be approximated by t'^'(A) oc (1/A)p. At the bot- 
tom of minimum light, p is approximately 0, and changes 
to 1, as the star returns to maximum light. By the 
least square fit for the extinction curves shown in Fig. ^ 
we derive logT'^'(A) = 0.52 — 0.47 log A(/im) when the 
V magnitude of the star is 10.61, while logT'^'(A) — 
-0.18-0.841ogA(/.tm) when the V magnitude is 7.03. The 
extinction at longer wavelengths is then estimated by ex- 
trapolation. Obviously this procedure is not optimal, but 
with no alternatives at hand for disentangling the contri- 
butions of the central star and the dust shell for the data 
studied here, we are forced to adopt this method. 

Assuming a blackbody for the flux of the central star, 
the flux density observed on the earth from the obscured 
central star is 



irBxiT^s) X exp(-T^'(A) - r^(A)), (5) 



where t^{X) is the optical depth of the dust shell along the 
radial direction. 

The flux density observed toward R CrB is, therefore, 
the sum of the flux from the central star obscured by a 
dust cloud and slightly dimmed by the optically thin dust 
shell (equation (||)), thermal emission from the optically 
thin dust shell (equation (||)), and the thermal emission of 
a newly formed (therefore presumably still optically thick) 
dust cloud given by 



/ci(A) 



^cl 

If 



^Ba(T,i) xexp(-r^(A)), 



(6) 



where Rd and Td are the radius and the temperature of 
the newly formed dust cloud, respectively. 
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The intensity distribution for the central star and the 
dust shell can be written as 

Ix{b) = BA(T.)exp(-r^'(A) - t%X)) x circ(6/i?0 



2mdKA / n{r)Bx{T{r))dz 



BA(Tci)exp(-T^(A)) X circ(&/i?ci), 



(7) 



where b is the impact parameter and z = ^/r"^ — 6^. The 
function circ(6/i?^) takes a value of 1 for |6| < i?* and 
elsewhere. The lower limit of the integration is Zmin — 

for 6 > Tin, while Zmin = V ~ ^^'^ ^ <^ ^in- The 
visibility is calculated by taking the modulus of the Fourier 
transform of the intensity distribution. 

5. Simultaneous fit of observed visibilities and 
SEDs 

5.1. Model fitting without thermal emission from a 
dust cloud 

We first try to fit the observed SEDs and visibilities us- 
ing models without thermal emission from a dust cloud, 
namely, neglecting the term given by equation (^). We 
adopt 7 = 2, appropriate for a constant mass loss and ex- 
pansion velocity. The sensitiveness of the near- and mid- 
infrared parts of emergent SEDs to Tin allows us to deter- 
mine it by fitting the observed SEDs. The visibility func- 
tions, which directly refiect the spatial extent of the shell, 
enable us to examine the validity of the models described 
above. 

Figs. ^ and ^ show the fit of the observed SEDs and 
visibility functions at the two epochs of our speckle obser- 
vations. The SEDs and visibilities are calculated with the 
opacities of amorphous carbon obtained by Busso letti et 
al. ( |1987D (AC2 sample). Rouleau & Martin (|l99l| ) (ACl 
sample), and Colangeli et al. ( 1995| ) (AGAR sample). The 
corresponding temperature distributions shown in Fig. ^ 
are used in the calculations. Figs. |6^ and ||b demonstrate 
that the SED and visibility for the 1996 data are well re- 
produced with Tin = 76 - 90 i?* and Tin = 900 - 950 K. 
For a given opacity data set, the uncertainty of Tin is esti- 
mated to be ±100 K, translating to an uncertainty of Hn 
of ±20 i?*. Taking the average of the Tin and ^n derived 
with three different opacities and adding the uncertainties 
resulting from the fitting, we derive Tin = 920 ± 103 K 
and Tin = 82 ± 23 R-i,. Fig. ^ illustrates the normalized 
intensity profile of the best fit model for the 1996 data. 
It consists of the central star and a ring-like structure 
characteristic of an optically thin shell. Note that the pre- 
dicted visibilities become plateau-like at spatial frequen- 
cies ^ 20 cycles/arcsec in Fig. ^jb. This plateau results 
from the unresolved central star. It should also be noted 
that the good match in the wavelength range shorter than 

1 fiui is simply due to the adoption of the empirical ex- 
tinction for a dust cloud, as described in Sect. |[ 

Figs. 0a and ^ reveal that the visibility observed 
at minimum light cannot be reproduced simultaneously 
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Fig. 6. Simultaneous fit of the SED and visibility of R 
CrB observed on 1996 October 1, using the models dis- 
cussed in Sect. 5.1. a The open squares, filled circles, and 
filled triangles represent the photometric data obtained 
by Shenavrin et al. ( 1979D , Shenavrin (unpublished obser- 
vation), and IRAS, respectively. The three curves repre- 
sent the SEDs predicted from models. CO represents the 
model with the data derived by Colangeli et al. (L995), 



RM with those derived by Rouleau & Martin (1991), and 
BU with those derived by Bussoletti et al. ( |1987| ). b The 
filled circles represent the observed visibility, while the 
three curves represent the visibilities predicted from the 
models, c Normalized intensity profile (A = 2.2 /im) of the 
BU model 



with the observed SED. The observed SED can be well 
fitted using the models with Tin = 700 - 800 K and ^n = 
116 - 172 i?*, but the visibilities predicted from these 
models are too low as compared with the observation. 
Regarding this discrepancy, we first examine the approx- 
imations adopted in our models. The line-blanketing ef- 
fect in the atmosphere of RCB stars is very prominent as 



Asplund et al. (1997) show, but the use of line-blanketed 
atmospheres instead of the blackbody would lower Tin by 
only ~ 50 K. The effect of the uncertainty of the effec- 
tive temperature is also minor. A decrease of the effective 
temperature by 500 K leads to a decrease of the dust tem- 
perature by ~ 60 K at ~ 100 i?*. Combining these effects. 
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Fig. 7. Simultaneous fit of the SED and visibility of R CrB 
observed on 1999 September 28, using the models with- 
out thermal emission from the newly formed dust cloud, 



as discussed in Sect. 5.1. a The filled circles represent the 
photometric data obtained one night after the speckle ob- 
servation. The IRAS data are represented by the filled 
triangles. The three curves represent the SEDs predicted 
from models. See also the legend to Fig. ^. b The filled 
circles represent the observed visibility, while the three 
curves represent the visibilities predicted from the models 



the dust temperature in the shell can be by ~ 100 K lower 
than those used in the fitting above. We tried to fit the 
observed SED and visibility using such a temperature dis- 
tribution, but it has turned out that the match to the 
observations is not much improved. 

One concern is the estimation of the /C-band flux of 
the central star obscured by the dust cloud. As described 
in Sect. ^ the extinction of the dust cloud in the JsT-band 
is derived by extrapolation from the region shortward of 
1 /im. However, the actual extinction curve of the dust 
cloud could become steeper longward of 1 /im and ap- 
proach the usual extinction curve of amorphous carbon, 
which is characterized by a spectral index of ~ 1.3 (Le 
Bertre |1997| ). We find that the observed SED and visibil- 
ity can be simultaneously fitted by adopting an extinc- 
tion curve for the obscuring dust cloud such as shown 
in Fig. ||c. The extinction curve bends at 1.4 /im and the 
spectral index changes from 0.47 to 1.3. Figs. |8^ and re- 
veal that the observed SED and visibility are fairly repro- 
duced. It should be stressed that the physical understand- 
ing of such an extinction curve is still unclear, and that 
the exact shape of the extinction curve cannot be uniquely 
determined by the fitting presented here. However, the un- 
usual properties of grains in the newly formed dust cloud 
may be responsible for the discrepancy found for the min- 
imum light data. 
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Fig. 8. Simultaneous fit of the SED and visibility of R 
CrB observed on 1999 September 28, using an extinction 
curve different from that extrapolated, a See the legend 
to Fig. 1^ for the reference of the symbols, b The filled 
circles represent the observed visibility, c The solid line 
represents an extinction curve for the dust cloud with a 
change of the spectral index at 1.4 /tm. The dashed line 
represents the extinction curve derived by extrapolation 
from the optical depth shortward of 1 /im, as discussed 
in Sect. ^. The filled squares represent the optical depths 
derived from the observations as shown in Fig. || 



5.2. Thermal emission from a newly formed dust cloud 

We show an alternative to interpret the data obtained at 
minimum light. We propose that the inconsistency found 
for the 1999 data set may be attributed to the thermal 
emission of a newly formed optically thick dust cloud in 
front of the star. The newly formed dust cloud is suppos- 
edly still rather close to the star and its angular size is 
not yet large enough to be resolved with the 6 m tele- 
scope. This assumption may be reasonable, because the 
1999 speckle observation was carried out in a decline when 
the star was obscured by AT^ = 4.6. The size and the tem- 
perature of such a dust cloud are by no means obvious, 
but the temperature (Td) should not exceed ~ 1500 K, 
otherwise the thermal emission of the cloud would be too 
prominent and would lead to a poorer match to the ob- 
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Fig. 9. Simultaneous fit of the SED and visibility of R CrB 
observed on 1999 September 28, using models with ther- 
mal emission from a newly formed dust cloud, as discussed 
in Sect. 5.2. a Three models with different data sets for 



the opacity of amorphous carbon are plotted. The radius 
and the temperature of a newly formed dust cloud are 
4.5 and 1200 K, respectively. CO, RM, and BU denote 
the opacities of amorphous carbon obtained by Colangeli 
et al. 
et al. 



m99q), Rouleau & Martin ( |199lD , and Bussoletti 
(1987), respectively. The contribution of the newly 
formed dust cloud is shown with the dash-dotted line for 
the BU model. See the legend to Fig. [t^ for the reference 
of the symbols, b The filled circles represent the observed 
visibility, c Normalized intensity profile (A = 2.1 /im) of 
the BU model. The model consists of the central star, a 
newly formed optically thick dust cloud with i?ci — 4.5 R^, 
(1.0 mas) and Tc\ = 1200 K, and an optically thin dust 
shell with ri„ = 180 R^, (41 mas) and Ti„ = 650 K 



served SED. We assume a circular cloud (projected onto 
the sky as seen from the star) with a radius of Rd be- 
tween 1 i?^ and 10 Ri,. As shown in Fig. ||, Rd denotes 
the radius of an individual cloud itself, not the distance be- 
tween the cloud and the central star. As 1 i?^, corresponds 
to 0.23 mas, even a cloud with a radius of 10 i?^ can well 
be treated as a point source for our speckle interferome- 
try. The thermal emission of the cloud is approximated by 
the blackbody radiation, as given by equation (||). Fig. ^ 
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Fig. 10. Simultaneous fit of the SED and visibility of 
R CrB observed on 1996 October 1, using models with 
thermal emission from a newly formed dust cloud out of 
the fine of sight, as discussed in Sect. 5^. A cloud with 
Rc\ = 4.5 Ri, and Td = 1200 K is placed out of the line 
of sight, 20 (5 mas) offset from the central star, a The 
three curves represent the SEDs predicted from the mod- 
els. CO, RM, and BU denote the opacities of amorphous 
carbon obtained by Colangeli et al. (1995), Rouleau & 
Martin (1991), and Bussoletti et al. ( 1987| ), respectively. 
The contribution of the newly formed dust cloud is shown 
with the dash-dotted line for the BU model. See the legend 
to Fig. ^ for the reference of the symbols, b The filled 
circles represent the observed visibility 



shows the normalized intensity profile of a model consist- 
ing of the central star, an optically thin dust shell, and a 
thermally emitting optically thick dust cloud, whose pa- 
rameters will be derived below. 

Figs. ^ and ||b show a simultaneous fit using models 
with thermal emission from a newly formed dust cloud. 
The observed SED and visibility are now fairly reproduced 
with models consisting of the central star, an optically 
thin dust shell with an inner radius of Tin — 150 - 193 i?* 
(average 172 R^) and Tin = 650 - 720 K (average 685 K), 
and an optically thick dust cloud with Rc\ = 4.5 i?* and 
Tci = 1200 K. The fit should be regarded as tentative. The 
uncertainties of Tin and Hn are estimated to be ±100 K 
and ±20 Ri,, respectively, while those of Td and Rd are 
about ±100 K and ±1 i?^, respectively. 

The difference of r-m derived for the 1996 and 1999 data 
may indicate a variation of the dust ejection frequency. 
If the star experiences dust ejection less frequently, the 
inner boundary is expected to become larger, and the IR 
excess is expected to decrease. In fact, as mentioned in 
Sect. |[ the L and M magnitudes decreased monotonically 
by ~ 1 mag from 1996 to 1999. This is consistent with the 
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trend found for rin derived for the two epochs, but further 
observations are required to confirm this correlation. 

Now we apply the thermally emitting cloud models to 
the data obtained near maximum light as well for the fol- 
lowing reason. Fig. ^ shows that R CrB started its final 
recovery from the 1995-96 deep minimum about 6 months 
before our speckle observation on 1996 October 1. If dust 
ejection occurs rather frequently, for example, every pul- 
sational cycle (40 - 50 days), we can expect that the star 
underwent dust ejection out of the line of sight during 
these 6 months. It is likely that newly formed dust clouds 
out of the line of sight existed on 1996 October 1, in ad- 
dition to the dust cloud in front of the star responsible 
for an obscuration of AV^ ^ 1. The latter cloud may al- 
ready be regarded as part of the optically thin dust shell, 
because = 0.34 is derived for this cloud by the extrap- 
olation described in Sect. ^ Therefore, for near maximum 
light, we consider a model with the central star dimmed 
by Ay ~ 1, thermal emission from an optically thin dust 
shell, and that from a newly formed optically thick dust 
cloud out of the line of sight. Fig. shows the SED and 
visibility calculated with an optically thick dust cloud with 
a radius of Rd — 4.5 i?* and Tci = 1200 K, in addition 
to the dimmed central star and the optically thin dust 
shell, whose parameters are the same as shown in Fig. ^. 
The cloud is placed out of the line of sight, 20 i?^ (5 mas) 
offset from the central star in the plane of the sky. Since 
the central star has already almost regained its brightness 
near maximum light, the effect of the newly formed dust 
cloud is minor on the SED and visibility. Therefore, the 
models including thermal emission from a newly formed 
dust cloud can also provide good matches to the SED and 
visibility observed near maximum light. 

6. Concluding Remarks 

Our 75 mas resolution speckle interferometric observations 
with the SAO 6 m telescope have spatially resolved the 
dust shell around R CrB for the first time. Neither the 
observation near maximum light nor at minimum light 
shows any clear deviation from circular symmetry. 

In order to derive the size of the dust shell, we first 
considered models consisting of the central star and an 
optically thin dust shell, neglecting the thermal emission 
of a newly formed dust cloud. Simultaneous fits of the 
models to the observed SED and visibility have demon- 
strated that a model with the central star and an optically 
thin dust shell with density proportional to r^^ seems to 
be appropriate for R CrB near maximum light. The inner 
boundary is found to be 82 R^, (19 mas) ± 23 R^, with a 
temperature of 920 ± 103 K. 

This simple picture fails to simultaneously reproduce 
the SED and visibility observed at minimum light, which 
has led us to investigate models with thermal emission 
from a newly formed optically thick dust cloud whose an- 
gular size is not yet large enough to be spatially resolvable. 
The SED and visibility obtained at minimum light were 
shown to be well fitted with such models. The presence of 



a newly formed dust cloud as hot as 1200 K with a radius 
of 4 - 5 Ri, is inferred, together with an optically thin dust 
shell with rjn ~ 170 i?* and Tin ^ 690 K. Furthermore, 
the SED and visibility obtained near maximum light were 
shown to be fitted also using a model with a newly formed 
dust cloud out of the line of sight. However, we have also 
discussed that the discrepancy found for the minimum 
light data may be attributed to the unusual extinction 
curve of the obscuring dust cloud. Observations during 
the very bottom of a deep minimum, when the contribu- 
tion of the central star is truly negligible, are crucial for 
investigating the dust shell and dust clouds. 
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